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Abstract. Eigenvectors of tensors, as studied recently in numerical multilinear algebra, 
correspond to fixed points of self-maps of a projective space. We determine the number of 
eigenvectors and eigenvalues of a generic tensor, and we show that the number of normal- 
ized eigenvalues of a symmetric tensor is always finite. We also examine the characteristic 
polynomial and how its coefficients are related to discriminants and resultants. 



1. Introduction 

In the current numerical analysis literature, considerable interest has arisen in extending 
concepts that are familiar from linear algebra to the setting of multilinear algebra. One such 
familiar concept is that of an eigenvalue of a square matrix. Several authors have explored 
definitions of eigenvalues and eigenvectors for higher- dimensional tensors, and they have 
argued that these notions are useful for wide range of applications [8]. Our approach in this 
paper is based on the definition of E-eigenvalues of tensors introduced by Liqun Qi in [9| [XI]. 
Throughout this paper, eigenvalue will mean E-eigenvalue, as defined in Definition 11.11 

We fix two positive integers m and n, and we consider order- m tensors A = (a^...^) of 
format n x n x • • ■ x n with entries in the field of complex numbers C. 

Definition 1.1. Let x be in C n and A a tensor as above. Adopting the notation introduced 
in [9j [TTJ, we define Ax m ~ l to be the vector in C n whose j-th coordinate is the scalar 



12=1 «m = l 

„m— 1 



If A is a complex number and x G C ra a non-zero vector such that Ax m = Ax, then A is 
an eigenvalue of A and x is an eigenvector of A. We will refer to the pair of A and x as 
an eigenpair. Two eigenpairs (A, x) and (A', x') of the same tensor A are considered to be 
equivalent if there exists a complex number t ^ such that t m ~ 2 X = A' and tx = x'. 

If m = 2 then Ax 1 is the ordinary matrix-vector product, and Definition 11.11 recovers 
the familiar eigenvalues and eigenvectors of a square matrix A. In that case, our notion of 
equivalence amounts to rescaling the eigenvector, but the eigenvalue is uniquely determined. 

For m > 3, Qi [TTJ normalizes the eigenvectors x of A by additionally requiring x ■ x = 1. 
When x ■ x is non-zero, this has the effect of choosing two distinguished representatives, 
related by A' = (— l) m A and x' = —x, from each equivalence class. In particular, when m 
is even, the eigenvalue is uniquely determined, and when m is odd, it is determined up to 
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sign. However, since equivalence classes with x-x = are not allowed by Qi's normalization, 
his definition does not strictly generalize the classical eigenvalues of a matrix. We will call 
eigenvalues A with an eigenvector satisfying x-x = 1 normalized eigenvalues of the tensor A. 

In Section 6 of [11], Qi considers an alternative normalization by requiring x-x = 1, where 
x is the complex conjugate. This reduces the equivalence classes from two real dimensions 
to one real dimension. One can still get an equivalent eigenpair (t m ~ 2 \,tx) for any complex 
number t with unit modulus. Yet another normalization, based on the p-norm over the real 
numbers R, was introduced by Lek-Heng Lim in his variational approach [8]. 

For most of this paper, we prefer not to choose any normalization whatsoever. Instead, 
we depend on the notion of equivalence in Definition 11.11 in order to have a finite number 
of equivalence classes of eigenpairs in the generic case. This equivalence is a generalization 
of the usual ambiguity of eigenvectors of an nxn-matrix A, which at best, are only unique 
up to scaling. The following theorem generalizes, from m = 2 to m > 3, the familiar linear 
algebra fact that an nxn- matrix A has precisely n eigenvalues over the complex numbers. 

Theorem 1.2. If a tensor A has finitely many equivalence classes of eigenpairs over C 
then their number, counted with multiplicity, is equal to ((m — l) n — l)/(m — 2). // the 
entries of A are sufficiently generic, then all multiplicities are equal to 1, so there are exactly 
((m — l) n — l)/(m — 2) equivalence classes of eigenpairs. 

For the case when the tensor order m is even, the above formula was derived in [9, Theo- 
rem 3.4] by means of a detailed analysis of the Macaulay matrix for the multivariate resultant. 
For arbitrary m, it was stated as a conjecture in line 2 on page 1228 of (9]. We here present 
a short proof of this conjecture that is based on techniques from toric geometry [21 H]. 

The rest of this paper is organized as follows. In Section [21 we prove Theorem 11.21 In 
Sectional we investigate the characteristic polynomial. Tensors with vanishing characteristic 
polynomial are interpreted as singular tensors. In Section HI we relate eigenvalues of tensors 
to dynamics on projective space. Finally, in Section [SI we specialize to the case of symmetric 
tensors. We show that, in that case, the set of normalized eigenvalues is always finite. 

2. Intersections in a Weighted Projective Space 

We shall formulate the problem of computing the eigenvalues and eigenvectors of the 
tensor A as an intersection problem in the n-dimensional weighted projective space 

X = P(l,l,...,l,m-2). 

The textbook definition of X can be found, for example, in [21 §2.0] and [U page 35]. 
Points in X are represented by vectors of complex numbers (ui : ■ ■ ■ : u n : A), not all zero, 
modulo the rescaling (tui : • • • : tu n : t m ~ 2 X) for any non-zero complex number t. The 
corresponding algebraic representation of our weighted projective space is X = Proj(-R), 
where R = C[xi, . . . ,x n , A] is the polynomial ring with x±, . . . ,x n having degree 1 and A 
having degree m — 2. The following proof uses basic toric intersection theory as in [H Ch. 5]. 
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Proof of Theorem \1.2 . For m = 2, the expression ((m — 1)™ — l)/(m — 2) simplifies to n, 
which is the number of eigenvalues of an ordinary nxn- matrix. Hence we shall now assume 
that m > 3. For a fixed tensor A, the n equations determined by Ax" 1 " 1 = Xx correspond 
to n homogeneous polynomials of degree m — 1 in our graded polynomial ring R. 

Since R is generated in degree m — 2, the line bundle Ox{m — 2) is very ample. The 
corresponding lattice polytope A is an n- dimensional simplex with vertices at (m — 2)ej for 
1 < % < n, and e n+ i, where the e, are the basis vectors in M n+1 . The affine hull of A is the 
hyperplane x\ + • • • + x n + (m— 2)A = m— 2. The normalized volume of this simplex equals 



(2.1) Vol(A) = (m-2 



,71-1 



The lattice polytope A is smooth, except at the vertex e n+ i, where it is simplicial with index 
m — 2. Therefore, the projective toric variety X is simplicial, with precisely one isolated 
singular point corresponding to the vertex e n+1 . By [U p. 100], the variety X has a rational 
Chow ring A*(X)q, which we can use to compute intersection numbers of divisors on X. 

Our system of equations Ax m ~ l = Xx consists of n polynomials of degree m — 1 in R. Let 
D be the divisor class corresponding to Ox{m — 1), and let H be the very ample divisor class 
corresponding to Ox(m — 2). The volume formula (12.11) is equivalent to H n = (m — 2) n ~ 1 in 
A*(X)q, and we compute the self-intersection number of D as the following rational number: 

D n = -ifY = ( r ^l) n .(m-2r- 1 - (m ~ 1)n 



m — 2 / \m — 2/ m — 2 

From this count we must remove the trivial solution {x = 0} of Ax 771 " 1 = Xx. That 
solution corresponds to the singular point e n+ \ on X. Since that point has index m — 2, the 
trivial solution counts for l/(m — 2) in the intersection computation, as shown in [H p. 100]. 
Therefore the number of non-trivial solutions in X is equal to 

(2.2) = 

v ; m-2 m-2 

Therefore, when the tensor A admits only finitely many equivalence classes of eigenpairs, 
then their number, counted with multiplicities, coincides with the positive integer in ( 12. 21) . 

In Example !2.2l below we exhibit a tensor A which attains the upper bound (12.21) . For that 
A, each solution (x, X) has multiplicity 1. It follows that the same holds for generic A. □ 

Remark 2.1. An alternative presentation of our proof is to perform the substitution A = 
\m-2 j n equations Ax 171 " 1 = Xx. This makes the system of equations homogeneous of 
degree m — 1. Bezout's theorem says that there are generically (m — l) n solutions in the 
projective space P n . If we remove the trivial solution, this leaves (m — l) n — 1 solutions. 
Assuming that none of these have A = 0, the orbits formed by multiplying A by e 2m /( m - 2 ) 
each yield the same value of A = X m ~ 2 and x. Thus, there are ((m — l) n — l)/(m — 2) classes. 

The delicate point in such a proof would be to argue that the solution to Ax 771 ^ 1 = 
X m ~ 2 x has multiplicity m — 2 even when A = 0. In effect, toric geometry conveniently does 
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the bookkeeping in the correspondence between solutions to Ax m 1 = Ax and solutions to 
Ax m-i = \ m ~ 2 x. 

Example 2.2. Let A be the diagonal tensor of order m and size n defined by setting Aa___i = 1 
and all other entries zero. An eigenpair (X,x) is a solution to the equations 

(2.3) x™- 1 = \xi for 1 < i < n. 

All non-trivial solutions in X = P(l, . . . , 1, m — 2) satisfy A ^ 0. By rescaling, we can assume 
that A = 1. Fix the root of unity ( = e 2 W( m - 2 ) ^ anc j g _ . . . ; 777. — 3, *}. For any string 
a in S n other than the all *s string, we define x, = ( ai if <jj is an integer and x, = if o"j = *. 
This defines (m — l) n — 1 eigenpairs. However, some of these are equivalent. Incrementing 
each integer in our string by 1 modulo m — 2 corresponds to multiplying our eigenvector 
by (. Thus, we have defined ((m — l) n — l)/(m — 2) equivalence classes of eigenvalues and 
eigenvectors. These are all equivalence classes of solutions to (I2.3p . 

More generally, suppose that A is a diagonal tensor with Au___i equal to some non-zero 
complex number a^. Then the eigenpairs are similarly given by A = 1 and Xj = aj m 2 ( l Ci 
or Xj = 0, as above, where aj m ~' 1 is a fixed root of a t . In particular, for generic aj all 
((m — l) n — l)/(m — 2) eigenpairs will have distinct normalized eigenvalues. □ 

From Theorem ll.2[ we get the following result guaranteeing the existence of real eigenpairs. 

Corollary 2.3. If A has real entries and either m or n is odd, then A has a real eigenpair. 

Proof. When either m or n is odd, then one can check that the integer ((to— l) n — l)/(m — 2) 
in Theorem 11.21 is odd. This implies that A has a real eigenpair by [3j Corollary 13.2]. □ 

Corollary 12.31 is sharp, in the sense that there exist real tensors with no real eigenpairs 
whenever both m and n are even. We illustrate this in the following example. 

Example 2.4. Let m be even, n = 2, and A the 2x • • • x2 tensor which is zero except for 
the entries ai 2 ...2 = 1 and a 2 i...i = — 1. The eigenpairs of A are the solutions to the equations: 

x™- 1 = \ Xl 
-x™- 1 = Ax 2 . 

Eliminating A, we obtain x™ + x?f = 0, which has no non-zero real solutions for even m. 

For n any even integer, let B be the tensor whose n/2 diagonal 2 x • • • x 2 blocks are the 
tensor A above, and which is zero elsewhere. A non-trivial eigenpair must be an eigenpair for 
at least one of the blocks, and therefore cannot be real. Thus, B has no real eigenpairs. □ 

3. Characteristic Polynomial and Singular Tensors 

The characteristic polynomial (ft a(A) of a generic tensor A was defined as follows in P| ITTj. 
Consider the univariate polynomial in A that arises by eliminating the unknowns 
from the system of equations Ax m ~ x = Ax and x • x = 1. If m is even then this polynomial 
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equals 0a(A). If m is odd then this polynomial has the form 0a(A 2 ), i.e. the characteristic 
polynomial evaluated at A 2 . With these definitions, Theorem 11.21 implies the following: 

Corollary 3.1. For a generic tensor A, the characteristic polynomial <pA (A) is irreducible 
and has degree ((m — l) n — l)/(m — 2). Hence this is the number of normalized eigenvalues. 

For any particular tensor A, the characteristic polynomial 0a(A) is obtained by specializing 
the entries in the coefficients of the generic characteristic polynomial. Ni et al. [9] expressed 
^(A) as a Macaulay resultant, which implies a formula as a ratio of determinants. For the 
present work, we used Grobner-based software to compute the characteristic polynomials 
of various tensors. It is tempting to surmise that all zeros of the characteristic polynomial 
4>a(X) are normalized eigenvalues of the tensor A. This statement is almost true, but not 
quite. There is some subtle fine print, to be illustrated by Example 13.21 below. 

Qi [Tlj Question 1] asked whether the set of normalized eigenvalues of a tensor is either 
finite or all of C. We answer this question by showing a tensor where neither of these 
alternatives holds: 

Example 3.2. Consider the complex 2x2x2 tensor A whose nonzero entries are 

sui = a-221 = 1 and a U2 = a 222 = i = v / -T- 

We claim that any complex number other than is a normalized eigenvalue of A, but is 
not a normalized eigenvalue. The equations for an eigenvalue and eigenvector of A are 

x\ + ix\X 2 = Xxi and X\X 2 + %x\ = Xx 2 . 

For any A ^ we obtain a matching eigenvector that also satisfies x ■ x = 1 by taking 

V 2A ' 2*A )' 

Hence A is a normalized eigenvalue. However, if A = 0, then an eigenvector must satisfy 

x\ + ix\x 2 = and xix 2 + ix\ = 0. 

These imply that x ■ x = x\ + x\ is zero, so A = cannot be a normalized eigenvalue. □ 

However, we have the following weaker statement: 

Proposition 3.3. The set of normalized eigenvalues of a tensor is either finite or it consists 
of all complex numbers in the complement of a finite set. 

Proof. The set S(A) of normalized eigenvalues A of the tensor A is defined by the condition 

3 x e C n : Ax 171 ' 1 = Xx and x ■ x = 1. 

Hence S(A) is the image of an algebraic variety in C n+1 under the projection (x, A) 1— > X. 
Chevalley's Theorem states that the image of an algebraic variety under a polynomial map 
is constructible, that is, defined by a Boolean combination of polynomial equations and 
inequations. We conclude that the set S(A) of normalized eigenvalues is a constructible 
subset of C. This means that S(A) is either a finite set or the complement of a finite set. □ 
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The relationship between the normalized eigenvalues and the characteristic polynomial is 
summarized in the following proposition. 

Proposition 3.4. For a tensor A, each of the following conditions implies the next: 

(1) The set £(A) of all normalized eigenvalues consists of all complex numbers 

(2) The set £(A) is infinite. 

(3) The characteristic polynomial </>a(A) vanishes identically. 

Proof. Clearly, (1) implies (2). By the projection argument in the proof above, the zero set 
in C of the characteristic polynomial 0a(A) contains the set £{A). Hence (2) implies (3). □ 

From Example 13.21 we see that (2) does not necessarily imply (1), and in Example 13.51 we 
shall see that (3) does not imply (2). Qi defines a singular tensor to be one for which the first 
statement of Proposition 13.41 holds. However, we suggest that the last condition is a better 
definition: a singular tensor is a tensor such that 0a(A) vanishes identically. This definition 
has the advantage that the limit of singular tensors is again singular. In particular, the set 
of all singular tensors is a closed subvariety in the n m -dimensional tensor space C nx '" xn . 
Its defining polynomial equations are the coefficients of the characteristic polynomial <pA (A) 
where A is the tensor whose entries a^...^ are indeterminates. 

Example 3.5. Let m = 3 and n = 2. Here A = (aijk) is a general tensor of format 2x2x2. 
The characteristic polynomial <$a is obtained by eliminating X\ and X2 from the ideal 

{aiiixl + {aii2+ai2i)xix 2 + ai22x\-\xi , a2iix\ + (a2i2 + ai2i)xiX2+a222xl-\x2 , x\+x\-l ). 

We find that 4>a has degree 3, as predicted by Theorem 11.21 Namely, the elimination yields 

0a(A 2 ) = C 2 A 6 + C 4 A 4 + C 6 A 2 + C 8 , 

where Cj is a certain homogeneous polynomial of degree i in the eight unknowns o^-fc. The 
set of singular 2 x 2 x 2-tensors is the variety in the projective space P 7 = P(C 2x2x2 ) given by 

(3.1) (C 2 ,C 4 ,C 6 ,C 8 ) c C[o m 

This is an irreducible variety of codimension 2 and degree 4, but the ideal (13. ip is not prime. 
The constant coefficient Cg is the square of a quartic. That quartic is the Sylvester resultant 

( dm a 112 + a 122 \ 

dill 0-112 + 0-121 G-122 

G-211 O212 + O22I 0-222 

\ 0-211 0-212 + O22I O222 / 

The leading coefficient of the characteristic polynomial is a sum of squares 

C2 = (— 0-111 + Ol22 + O212 + 022l) 2 + (O112 + O121 + O211 — 0222) 2 

This indicates that among singular 2 x 2 x 2-tensors those with real entries are scarce. Indeed, 
the real variety of (13.11) is the union of two linear spaces of codimension 4, with defining ideal 

(Ol22, 0-211) O112 + O121 — 0222) O212 +O22I — Oin) fl (a m — 0-122, O211 — 0222) ^112 + 0121) a 212 +022l)- 



Res x (Ac 2 ) = det 
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This explains why the singular tensor in Example 13.21 had to have a non-real coordinate. 

We now look more closely at the real singular tensors defined by the second ideal in this 
intersection. These are tensors A for which Ax 2 = (am(x 2 + y 2 ) , a 2 ii(x 2 + y 2 ))- It is easy 
to see that, so long as am and 0211 are not both zero, the only normalized eigenvector is 

(am 0211 
\Am + a ln V a in + a lu 

In particular, the number of eigenvalues of such a tensor must be finite. This example shows 
that (3) does not imply (2) in Proposition 13.41 □ 

The reader will not have failed to notice that the notion of "singular" used here (and 
in [11]) is more restrictive than the one familiar from the classical case m — 2. Indeed, a 
matrix is singular if it has A = as an eigenvalue, or, equivalently, if the constant term 
of the characteristic polynomial vanishes. That constant term is a power of the resultant 
ReSa;(y4x m_1 ), and its vanishing means that the homogeneous equations Ax 171 ' 1 = have a 
non-trivial solution x G P" -1 . This holds when the tensor A is singular but not conversely. 

Yet another possible notion of singularity for a tensor A arises from its hyperdeterminant 
Det(A), as defined in [5]. For example, the hyperdeterminant of a 2x2x2-tensor equals 

Det^) = a 122 a 211 a 121 a 212 a 112 a 221 a lll°222 

— 2a 1 2ia.i22 c !'2ii a 2i2 — 2a 112 ai22 a 2ii a 22i — 2a 112 ai2ia2i2022i — 2a 111 a 12 2a-2iifl222 
— 2aniai2ia2i20222 — 2ama112a.221a.222 + 4amai22a2i2a22i + 4an2ai2ia2na222- 

The hyperdeterminant vanishes if the hypersurface defined by the multilinear form associated 
with A has a singular point in (p n_1 ) m , This property is unrelated to the coefficients of the 
characteristic polynomial 0a(A). In particular, Det(v4) 7^ Res^Ac" 1 " 1 ) . 

For an example take the 2x2x2-tensor A all of whose entries are = 1 is not singular 
but Det(v4) = 0. On the other hand, the following tensor A is singular but has Det(v4) = —1: 

am = —1, an2 = 0, am = 0, ai22 = —1, 0211 = 1, 0212 = —1, 0221 = 0, 0222 = 1 — 
This highlights the distinction between our setting here and that in [HI Proposition 2]. 

4. Dynamics on Projective Space 

The purpose of this short section is to point out a connection to dynamical systems. 
Dynamics on projective space is a well-established field of mathematics [TJ [7J. We believe 
that the interpretation of eigenpairs of tensors in terms of fixed points on P n_1 could be of 
interest to applied mathematicians as a new tool for modeling and numerical computations. 

We consider the map if) a defined by the formula iPa{x) = Ax™ 1 ' 1 . This is a rational map 
from complex projective space P n_1 to itself. The fixed points of the map if) a '■ P n_1 — - » P™ -1 
are exactly the eigenvectors of the tensor A with non-zero eigenvalue, and the base locus 
of if) a is the set of eigenvectors with eigenvalue zero. In particular, the map ipA is defined 
everywhere if and only if is not an eigenvalue of A. 



which has eigenvalue A = \ / a m + a\ xl 
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Note that every such rational map arises from some tensor A, but the tensor is not unique. 
Indeed, A has n m entries while the map is determined by n polynomials, which have only 
n("] i ™j 2 ) distinct coefficients. For instance, in Example I3.5[ with m = 3,n = 2, the eight 
entries of the tensor translate into six distinct coefficients of the two binary quadrics that 
specify the self-map of the projective line ipA'- P 1 -—* P 1 - 

It is instructive to revisit the classical case m — 2, where ipA'- P n_1 — •* P"-^ 1 is a linear 
map. The condition that every eigenvalue of the matrix A is zero is equivalent to saying 
that A is nilpotent, that is, some matrix power of A is zero. Geometrically, this means that 
some iterate of the rational map ipA is defined nowhere in projective space P n_1 . We use the 
same definition for tensors: A is nilpotent if some iterate of i/ja is nowhere defined. 

Proposition 4.1. If the tensor A is nilpotent then is the only eigenvalue of A. The 
converse is not true: there exist tensors with only eigenvalue that are not nilpotent. 

Proof. Suppose A 7^ is an eigenvalue and x G C n \{0} a corresponding eigenvector. Then x 
represents a point in P n_1 that is fixed by ipA- Hence it is fixed by every iterate ip^ °f V'a- 
In particular, if)^ is defined at (an open neighborhood) of x G P n_1 , and A is not nilpotent. 

Let A be the 2x2x2-tensor with a m = a 2 n — ^212 — 1 and the other five entries zero. 
The eigenpairs of A are the solutions to x\ = Xxi and x\ + x±x 2 = Xx 2 . Up to equivalence, 
the only eigenpair is x = (0, 1) and A = 0. However, the self-map ipA on P 1 is dominant. 
To see this, note that ipA acts by translation on the affine line A 1 = {xi 7^ 0} because 
(x 2 : x\ + X1X2) = (xi : xi + X2). All iterates of if) a are defined on A 1 , i.e. there are no base 
points with x\ 7^ 0, and hence A is not nilpotent. □ 

The example in the previous proof works because the two binary quadrics in Ax m ~ l have 
X\ as a common factor. Indeed, whenever n = 2, an eigenvector x has eigenvalue zero if and 
only if £ is a solution to a common linear factor of the two binary forms of Ax" 1-1 . 

However, for n > 3, this is no longer true. Work in dynamics by Ivashkovic [TJ Theorem 1] 
implies that, for n — 3, one can construct tensors A such that zero is the only eigenvalue, 
the polynomials in Ax m ~ x have no common factors, but A is not nilpotent. 

Example 4.2. This example is taken from [7J Example 4.1]. Let m = n = 3 and take A to 
be any tensor whose corresponding map is the Cremona transformation 

ip A - P 2 — » P 2 : (xi,x 2 ,x 3 ) t-> (xiX 2 ,XiX 3 ,2x 2 X3). 

This map has no fixed points, but it is not nilpotent. The base locus of ipA consists of the 
three points (1, 0, 0), (0, 1, 0), and (0, 0, 1), and, up to scaling, these are the only eigenvectors 
of A, all with eigenvalue 0. □ 

5. Symmetric Tensors 

Of particular interest in numerical multilinear algebra is the situation when the tensor A 
is symmetric and has real entries. Here A being symmetric means that the entries a^...^ are 
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invariant under permuting the n indices i\, 12, ■ ■ ■ , i n - Each symmetric nx ■ ■ ■ xn tensor A of 
order m corresponds to a unique homogeneous polynomial f(x) of degree m in n unknowns. 

The symmetric case is of interest because a real polynomial f(x) of even degree m is posi- 
tive semidefmite if and only if every real eigenpair of the corresponding symmetric tensor A 
has non-negative eigenvalue [TUJ Theorem 5(a)]. This is illustrated in Example 15.81 

In the notation of [TO], the relation between the tensor and the polynomial is written as 

(5.1) Ax m = mf(x) and Ax 171 ' 1 = V/(x), 
where Ax m is defined to be 

n n n 

^ 1 ^ ] ' ' ' ^ 1 a ii-..i m x ii x h ' ' ' x im = x ' Ax 

11 = 142=1 im=l 

The first equation in (15. ip follows from the second because x ■ V/(x) = mf(x). Note that 
the second equation in (15. ip says that the coordinates of the gradient of f(x) are precisely 
the entries of Ax m ~ x . The gradient V/(x) vanishes at a point x in P n_1 if and only if x is a 
singular point of the hypersurface in P n_1 defined by the polynomial f(x). This implies: 

Corollary 5.1. The singular points of the projective hypersurface {x £ p™- 1 : /(ar) = 0} are 
precisely the eigenvectors of the corresponding symmetric tensor A which have eigenvalue 0. 

The other eigenvectors of A can also be characterized in terms of the polynomial f(x). 

Proposition 5.2. Fix a non-zero A and suppose m > 3. Then x £ C n is a normalized 
eigenvector with eigenvalue A if and only if x is a singular point of the ajfine hypersurface 
defined by the polynomial 

(5.2) /W _^. X _(I_1) A . 

Proof. The gradient of (15. 2p is V/ — Ax = Ax m ~ x — Ax, so every singular point x is an 
eigenvector with eigenvalue A. Furthermore, if we substitute f{x) = —x ■ V/ = ^x ■ x into 
(15. 2p . then we obtain x-x = 1. This argument is reversible: if x is a normalized eigenvector of 
A then x ■ x = 1 and V/(x) = Ax, and this implies that f 1 5 . 2 1) and its derivatives vanish. □ 

Corollary 5.3. The characteristic polynomial </>a(A) is a factor of the discriminant of 

Here we mean the classical multivariate discriminant [5] of an inhomogeneous polynomial 
of degree m in n variables x evaluated at (15. 2p . where A is regarded as a parameter. Besides 
the characteristic polynomial 0a(A), this discriminant may contain other irreducible factors. 

Example 5.4 (Discriminantal representation of the characteristic polynomial of a symmetric 
tensor). If n = 2 and m = 3 then the discriminant at bivariate cubic in (15. 2p equals A 4 -0 J 4(A). 
If n = 2 and m = 4 then we evaluate the discriminant of the ternary quartic using Sylvester's 
formula jS] §3.4.D]. The output has the discriminant of binary quartic as an extraneous factor: 

Discriminant of (15. 2 p = (0^(A)) 2 • A 9 • Discriminant of f(x) 
It would be interesting to determine the analogous factorization for arbitrary m and n. □ 
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The subject of this paper is the number of normalized eigenvalues of a tensor. In Section 2 
we gave an upper bound for that number under the hypothesis that the number is finite. 
Remarkably, this hypothesis is not needed if we restrict our attention to symmetric tensors. 

Theorem 5.5. Every symmetric tensor A has at most ((m — 1)™ — l)/(m — 2) distinct 
normalized eigenvalues. This bound is attained for generic symmetric tensors A. 

Proof. It suffices to show that the number of normalized eigenvalues of every symmetric 
tensor A is finite. Recall from the proof of Theorem 11.21 that the set of eigenpairs is the 
intersection of n linearly equivalent divisors on a weighted projective space. Since these divi- 
sors are ample, each connected component of the set of eigenpairs contributes at least one to 
the intersection number. Therefore, the number of connected components of eigenpairs can 
be no more than ((m — l) n — l)/(m — 2). We conclude that the number of normalized eigen- 
values of A, if finite, must be bounded above by that quantity as well. Finally, Example 12.21 
shows that the bound is tight. 

We now prove that the number of normalized eigenvalues of a symmetric tensor A is finite. 
Let S be the affine hypersurface in C n defined by the equation x\ -\ — • + x\ = 1. We claim 
that a point x G S is an eigenvector of A if and only if x is a critical point of / restricted 
to S, in which case, the corresponding eigenvalue A equals —f(x). By definition, a point 
x G S is a critical point of f\s if and only if the gradient V(/|s) is zero at x. The latter 
condition is equivalent to the gradient V/ being a multiple of V(x\ + ■ ■ • + x 2 n — 1) = 2x. 
This is exactly the definition of an eigenvector. Finally, if x G S is a critical point of f\$, 
then mf(x) = x ■ V/(x) = Ax • x = A, and hence A = —f(x). 

Finally, to prove Theorem 15.51 we note that, by generic smoothness [HI Cor. III. 10.7], a 
polynomial function on a smooth variety has only finitely many critical values. Equivalently, 
Sard's theorem in differential geometry says that the set of critical values of a differentiable 
function has measure zero, so, by Proposition 13.31 that set must be finite. □ 

We note two subtleties about Theorem 15.51 First, it does not imply that the characteristic 
polynomial of every symmetric tensor is non-trivial. Second, the result is intrinsically tied 
to the normalization x ■ x = 1. We begin with an example of the first. 

Example 5.6. Let A be the symmetric 2x2x2 tensor with 

Qlll — — 2i , <2ii2 = 0\2\ = 0211 — 1 j a 122 = &212 = = , 0222 = 1- 

Then, up to equivalence, the only eigenvectors are (0, 1) with eigenvalue 1 and (1, i) with 
eigenvalue 0. Note that the second cannot be rescaled to be a normalized eigenvector, 
so the only normalized eigenvalue is 1. However, the characteristic polynomial of A is 
identically zero. The reason is that, for a small perturbation of A, the perturbation of the 
eigenvector (l,z) can take on any given normalized eigenvalue. □ 

No analogue of Theorem 15.51 is possible with the alternative normalization of requiring 
x ■ x — 1. In this case, each equivalence class yields infinitely many eigenvalues, which 
nonetheless have the same magnitude. However, the following example shows that the mag- 
nitudes of the eigenvalues with x ■ x — 1 may still be an infinite set. 
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Example 5.7. Let A be the symmetric 3x3x3 tensor whose non-zero entries are 

Gm = 2 and 0,122 = 0-212 = 0,221 = 0,133 = 0,313 = 0331 = ^ 
The eigenpairs of A are the solutions to the equations 

2xix 2 = Ax 2 , 
2x1X3 = AX3. 

For any aGC, the vector x = (1, za, a) is an eigenvector with eigenvalue A = 2. Rescaling, 
x/y/x ■ x is an eigenvector with unit length and eigenvalue 

2 

y/l + 2\a\' 

The magnitude of this eigenvalue can be any real number in the interval (0, 2 ]. Note that the 
family of eigenvectors above all have x-x = l,soA = 2is the only normalized eigenvalue. □ 

One application of eigenvalues of symmetric tensors is that these can be used to decide 
whether a polynomial / is positive semidefinite, i.e., whether f(x) > for all x G M. n . 

Example 5.8. The Motzkin polynomial f(x, y, z) = z & +x A y 2 +x 2 y A — 3x 2 y 2 z 2 is a well-known 
example of a positive semidefinite polynomial which cannot be written as a sum of squares. 
Let A be the corresponding 3x3x3x3x3x3-tensor. This tensor has 25 eigenvalues, counting 
multiplicities, six less than our upper bound of 31. Disregarding multiplicities, there are only 
four distinct eigenvalues. All four are real and they are equal to: (with multiplicity 14), 3/32 
(with multiplicity 8), 3/2 (with multiplicity 2), and 6 (with multiplicity 1). By [TOj Theorem 
5(a)], this confirms the fact that the Motzkin polynomial / is positive semidefinite. □ 
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